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In the first one of these two lectures, I give an introductory review of phase transitions in fi- 
nite temperature field theories. I highlight the differences between theories with global and local 
symmetries, and the similarities between cosmological phase transitions and phase transitions in su- 
perconductors. In the second one, I review the theory of defect formation in finite temperature phase 
transitions, emphasising the differences between systems with broken global and local symmetries. 
These same principles apply to relativistic and condensed matter systems alike. 



Lecture I: 

Phase Transitions in the Early Universe 



INTRODUCTION 

If one calculates the equilibrium properties of the Stan- 
dard Model at high temperatures, one finds several dif- 
ferent phases Below 100 MeV, the elementary 
particles behave more or less as at zero temperature. If 
we heat them up above this temperature, there is a tran- 
sition to a deconfined phase in which quarks and gluons 
are not confined into hadrons. At 100 GeV, there is an- 
other transition at which the weak and electromagnetic 
interactions unify. Extensions of the Standard Model, 
such as Grand Unified Theories (GUTs), typically pre- 
dict new phases at even higher temperatures. 

This phase structure is important in cosmology ill 
because the temperature of the universe was initially very 
high. It seems likely that had to be higher than 100 
GeV, because otherwise it is extremely difficult to ex- 
plain the baryon asymmetry of the universe — the fact 
that there is so much matter in the universe, but practi- 
cally no antimatter. As the universe expanded and cooled 
down to the current 2.7 K, it must have therefore under- 
gone at least the deconfinement and electroweak phase 
transitions mentioned above. Typical inflationary mod- 
els predict much higher temperatures, usually at least 
10^" GeV, and it is therefore likely that there were also 
other phase transitions. 

The cosmological consequences of these phase transi- 
tions depend sensitively on their properties. In many 
cases the consequences would be disastrous, and one can 
then rule out the theory that predicts the transition, 
or at least constrain its parameters. Examples of this 
are grand unified theories with strongly first-order phase 



transitions, as they would lead to a period of "old infla- 
tion" , and to a very inhomogeneous universe, or theories 
that predict massive domain walls or too many magnetic 
monopoles. In other cases, the effects of the phase transi- 
tion may be subtle, only barely observable, such as grav- 
itational waves, or even completely unobservable. Or, it 
may be that the consequences are absolutely crucial for 
our own existence, as in the theory of electroweak baryo- 
genesis which uses the electroweak phase transition 
to explain the baryon asymmetry. 

In this lecture, I will concentrate on those phase tran- 
sitions that can be understood as a spontaneous break- 
down of either a global or a local symmetry. This covers 
the electroweak and GUT transitions, and to a certain ex- 
tent also the deconfinement transition. I will also point 
out similarities to phase transitions in condensed matter 
systems such as superfiuids and superconductors. Un- 
like cosmological phase transitions, they can be studied 
experimentally and can therefore be used to test our the- 
oretical understanding of phase transitions in quantum 
field theories. 



DEFINITIONS 

The properties of a thermodynamic system in thermal 
equilibrium are given by the partition function Z (see 
Ref. Q), which can be generally written as 



Z = Ti p = Tr 



(1) 



'Lectures given at a COSLAB Workshop on "Cosmological Phase 
Transitions and Topological Defects" , 22-24 May 2003, Porto, Por- 
tugal. 



where H is the Hamiltonian of the system and /? = 1/kT 
is the inverse temperature. 

In many cases, it is more convenient to use its loga- 
rithm, which defines the free energy F = —TlnZ. This 
quantity has the property that its value is minimised in 
thermal equilibrium. 

Neither the partition function nor the free energy are 
observable quantities. However, any observable can be 
expressed as the derivative of F with respect to some 
external field. For instance, the expectation value of the 
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FIG. 1: The phase diagram of water (greatly simplified). 
The lines are first-order phase transitions. The transition 
line between water and vapour ends at a second-order critical 
point. 



operator X is defined as 

{X) = Z-^Tr Xp, 
and if we define an external field Jx by 



Z{Jx) = Tr e 



-t3H+JxX 



we find 



{X) = 



dJx 



InZ 



1 dF 
TdJ^ 



(2) 



(3) 



(4) 



A phase transition is defined as a point at which the 
free energy (or, equivalently, the partition function) is 
non-analytic. In general, it follows that at least some 
observables are non-analytic at the transition point and 
act as order parameters. 

The simplest type of phase transitions arc first-order 
transitions, in which some expectation value, say {X), is 
discontinuous. An everyday example of this is the tran- 
sition between water and ice. Let us assume that the 
transition takes place at temperature T^. Because of the 
discontinuity, the state of the system at depends on 
which side the critical temperature was approached from. 
Consequently, it is possible to have two phases coexist- 
ing at Tc, which is exactly what happens when we put 
ice into a glass of water. As the temperature is gradually 
changed through Tg, the discontinuous jump from one 
value of (X) to the other does not necessarily take place 
instantaneously at Tc- For instance, water can remain liq- 
uid in a supercooled state at temperatures below 0°C for 
a relatively long time, if it is pure enough and completely 
unperturbed. This happens in a meteorological condition 
known as freezing rain, in which the supercooled water 
freezes as soon as it hits the ground, causing hazardous 
driving conditions. More generally, this phenomenon is 
known as metastability. 



H 
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FIG. 2: The phase diagram of a typical ferromagnet. Because 
of the rotation symmetry, the line of first-order transitions is 
on the T axis. It ends in a second-order transition point at 
the Curie temperature Tc- The rotation symmetry is spon- 
taneously broken at T < To. 



Second-order phase transitions are those in which ev- 
ery expectation value is continuous, but some expecta- 
tion value {X) has a discontinuous derivative. This is 
a limiting case of a first-order transition in which the 
discontinuity is taken to zero, and therefore it typically 
requires fine tuning of some parameter. On a phase di- 
agram, second order transitions are typically end points 
of first order transition lines. One example of this is the 
critical point of water at which the first-order transition 
line between water and vapour ends. To reach this point, 
one has to tune both the temperature and the pressure to 
T = 647.096 K and p = 22.064 MPa. In contrast, reach- 
ing a first order transition point, one only has to tune 
the temperature provided that the pressure is below the 
critical value. 

There is, however, an exception to this rule. There 
are cases in which some symmetry forces the system to 
be on the first-order transition line. This happens, for 
instance, in a ferromagnet. If the temperature is below 
the Curie temperature, the magnetisation M becomes 
non-zero even if there is no external magnetic field H. 
Because the setting is rotation invariant, no direction of 
M is preferred. Therefore, the direction of M, which de- 
pends on how the external field was taken to zero, breaks 
the rotation invariance. This is known as spontaneous 
symmetry breaking. 

If we imagine changing the external field H gradually 
through zero, we see that M is discontinuous at _ff = 0. 
In other words, there is a first-order transition at H = 0. 
If we increase the temperature keeping H zero, we are 
therefore guaranteed to reach the second-order end point 
at the Curie temperature. 
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The above situation is very important in field theories, 
because they typically have a high degree of symmetry. 
As in the case of the ferromagnet, it often happens at low 
temperatures that the state of the system is not invariant 
under a symmetry of the theory and this typically leads 
to a second-order phase transition as the temperature is 
increased. 



GLOBAL SYMMETRY BREAKING IN FIELD 
THEORIES 

As a simple example of a field theory with spontaneous 
symmetry breaking, let us consider the Lagrangian 



1 



A 



(5) 



where is a real scalar field. The Lagrangian is sym- 
metric under (/) — > — 0, which correspond to the group 
Z2. 

The vacuum of the system is found by solving the 
Euler-Lagrange equations for a constant field, 

= 0, (6) 



which gives 



Ai 



= 0. 



(7) 



In addition to = 0, which is an unstable solution, we 
find two solutions, (\) ~ ±w, which are the possible vac- 
uum states of the system. Neither choice is invariant 
under the Z2 symmetry, and therefore the symmetry is 
spontaneously broken in the vacuum. 

At a non-zero temperature, the value of fluctuates 
from point to point. If these thermal fluctuations are 
small, we can write the field as the sum of the homo- 
geneous background contribution (^q and the fluctuation 



(x) = (/io + (50(x). 



(8) 



We assume that these fluctuations are Gaussian, which 
is equivalent to the one-loop approximation in perturba- 
tion theory, i.e., in the expansion in terms of Feynman 
diagrams. In many cases, this is a good approximation. 
The homogeneous mode (/)o is chosen in such a way that 
the mean value of the fluctuation in zero, ((5(/)) = 0. The 
variance, however, is non-zero, and has the value 



d^p 1 



1 



(27r)3 w e^" - 1 



(zero-temperature part), 
(9) 

where w = (p*^ + m^)^/'^. Here m is the mass of the fleld 
in the vacuum around which the expansion is done, and 
is given by the second derivative of the potential, 

= (10) 



The zero-temperature part is actually divergent, but 
we do not have to worry about that, because it is al- 
ready absorbed in the definition of the renormalised zero- 
temperature couplings. For T ^ m, Eq. ^ can be ex- 
panded as 



y2 

12 



mT 



(11) 



When the fluctuations 6(l> are taken into account, the 
homogeneous more (jjo feels an effective potential that is 
different from the "tree- level" potential in Eq. ((Sj. This 
is easiest to see by averaging Eq. (0 over fluctuations j3| , 



= \{{q^')-v^^)) 



(12) 



where we have used the fact that the fluctuations are 
symmetric and odd moments must therefore vanish. If we 
now use the leading term from Eq. l(TT)) . « T^/12, 

we flnd 



= 0. 



(13) 



This shows that the two non-zero minima only exist be- 
low a critical temperature Tc = 2v, and above that, the 
symmetry is restored. 

One can easily generalise this discussion to somewhat 
more complicated cases. If (pi is a multicomponent real 
scalar fleld with i S 1, . . . , A^, we can consider the La- 
grangian 



(14) 



where summation over the indices i is assumed. This La- 
grangian is invariant under SO(A^) rotations of the fleld 
(/),;. For two components, we can deflne a complex scalar 
fleld $ = {(1)1 + i(j>2), and we have 



A 

4 



(15) 



This Lagrangian is invariant under rotations of the com- 
plex phase angle of the scalar fleld i/i. , which means that 
it has a U(l) symmetry. Because of this symmetry, we 
can choose that the vacuum expectation value is real, and 
expand — (j)o + S(j)i and 02 = S(j)2- The variance of 
both 6(j)i and (502 are given by Eq. pi|l . The averaged 
Euler-Lagrange equation for the real part is 



= 



dC 

>501, 

= A((0?) + (020i)-„2(^^)) 

= X{<j)l + 3{S<Pl) + {6c^l)-v')<l)o. (16) 
Using Eq. (fTTT) . we flnd 



T2 



= 0, 



(17) 
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theory 


symmetry 


order parameter 


10^^ GeV 


GUT 


SU(5)^ SU(3)xSU(2)xU(l) 


<I>, 24 real components 


100 GeV 


electroweak 


SU(3)xU(l) 


<j), 2 complex components 


100 MeV 


QCD 


U(l)" 


qq composite 



TABLE I: The phase transitions predicted by particle physics models. 



showing that this time the critical temperature is Tc = 
-s/Sw, but otherwise the conclusions are unchanged. 

The restoration of broken symmetries at high temper- 
atures is important for the early universe, because the 
universe was initially extremely hot and because the the- 
ories of particle physics contain several symmetries that 
are spontaneously broken at zero temperature. As shown 
in Table these are associated with the hypothetical 
Grand Unified Theory (GUT), the electroweak theory 
and quantum chromodynamics (QCD). 

The QCD deconfinement phase transition takes place 
at around 100 MeV, and is the only one of the three tran- 
sitions that can be studied experimentally. The transi- 
tion cannot be fully understood within perturbation the- 
ory, and it can be characterised in two different ways, 
either as a transition from a confining to a deconfined 
phase, or as spontaneous breakdown of the chiral sym- 
metry. We will adopt the latter point of view. 

If the up and down quarks were massless, the left 
and right handed quarks would decouple. One could 
therefore carry out independent rotations between the 
up and down quarks of the given handedness. In the 
QCD vacuum, this symmetry is spontaneously broken 
by the quark condensate (gg), which shows how the left 
and right handed quarks are aligned in the flavour space. 
Just like in a ferromagnet, one would expect the symme- 
try to be restored at high temperatures, and this seems 
to imply a second-order phase transition. 

However, the quarks are not actually massless and 
therefore the chiral symmetry is not exact. The quark 
mass term plays the same role as the magnetic field H 
in the ferromagnet phase diagram in Fig. 13 and shifts 
the system a little bit away from the first-order transi- 
tion line. Because there is no symmetry in the first place, 
there is no second order transition either. A first-order 
transition is still possible if the phase diagram is more 
complicated than that in Fig. [3 

The GUT and electroweak phase transitions, on the 
other hand, involve exact symmetries, but local gauge 
symmetries instead of global ones. Therefore, the above 
discussion does not apply them either, as we shall see in 
the next section. 



point in order for the Lagrangian to remain unchanged. 
It is possible to modify the theory in such a way that it 
is invariant under position-dependent rotations as well. 
This is done by introducing a gauge field A^, and the 
local symmetries are therefore also known as gauge sym- 
metries. For the theory with a complex scalar in Ea. ll5|l . 
the modified Lagrangian is 

C - -\f^.F^'' + \d^^*D^^ - ^ {cj>*^ - v^f , (18) 

where we have used the covariant derivative I?^ = 9^ + 
igA^ and the field strength tensor = — d^A^j^. 

This model is essentially a relativistic version of the 
Ginzburg-Landau theory of superconductivity. Our re- 
sults will therefore apply to superconductors on a quali- 
tative level. The quantitative results differ, though. 

The gauge coupling constant g is simply the electric 
charge of the field 0. (For superconductors, g — — 2e, 
because a Cooper pair consists of two electrons.) The 
Lagrangian is invariant under local (position-dependent) 
rotations of the phase angle 



A„ 



A„ 



9 



(19) 



Again, we write = 0o + ^i'l where we now choose (po to 
be real and positive. The gauge field A^ does not have 
a homogeneous mode, and therefore we can treat is as a 
pure fluctuation. Again, we average the Euler-Lagrange 
equation for the real part 0i over fluctuations (5(/>i, (502 
and A„ , 



dC 



= A 



-v^^,))+g'{A^Af^^,) 



The gauge field contains three physical degrees of free- 
dom, and therefore we have 



{A.A'^) 



4 



47r 



(r = Opart). (21) 



GAUGE SYMMETRY BREAKING 



Crucially, however, the mass rn^ of the gauge field de- 
pends on 00 because of the Higgs mechanism, and is 
m~f — gcjjQ. Including this contribution, we find 



The Lagrangians discussed above had global symme- 
tries. One has to carry out the same rotation at every 



A 



4^T^^° 



3 4A/ 



= 0. (22) 
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FIG. 3: The critical temperature [solid line] in units of v 
at g = 0.5. The dashed lines show the metastability range. 



To interpret this results, we define the effective poten- 
tial as the tree-level potential that would give this same 
Euler-Lagrange equation. That is obtained by simply 
integrating the left hand side of the equation, 



An 



(23) 



The presence of the cubic term means that as long as 
the coefficient of the quadratic term is positive and small 
enough, there are two non-negative minima, one at (/iq — 
and the other at a non-zero value. The two minima are 
degenerate at a certain critical value of Tc, 



12Aw2 



4A + 3.g2- (3/27r2)(g6/A)' 



(24) 



This corresponds to a first-order transition line. Thus we 
see that due to thermal fluctuations, we have a first-order 
transition instead of a second-order one.^ The disconti- 
nuity in 00 is 



A0O 



1 

2^ A 



(25) 



It is also straightforward to calculate the range of tem- 
peratures in which a second, metastable minimum exists. 



12Aw2 
4A + 3g2 



4A + 352- (27/167r2)(g6/A) 



(26) 



The critical temperature Tc, together with the metasta- 
bility range is shown in Fig. O for g = 0.5. Fig. 0] shows 
the shape of the effective potential at Tc and at the ends 
of the metastability range for g — 0.5, A — 0.05. 



^ Note that this is the reason why Type I superconductors have a 
first-order transition. 



0.4 
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FIG. 4: The effective potential at Tc [solid line] for g = 0.5, 
A = 0.05. The dashed lines show the potential at the two 
ends of the metastability range. 



When deriving Eq. (|23|l . we have assumed that the fluc- 
tuations are Gaussian. It is important to check whether 
this approximation, which is equivalent to a one-loop cal- 
culation in terms of Feynman diagrams, is consistent. At 
higher orders in the loop expansion, one finds that each 
loop gives typically a contribution of the order g^T jm^. 
At Tc, the photon mass in the broken phase is 



(27) 



and therefore each loop order gives a contribution propor- 
tional toXj g^ . If A ^ ^2^ this is small and means that the 
contribution from higher loop diagrams can be ignored. 
One can therefore trust the calculation of the effective 
potential in Eq. I|23|l . On the other hand, if A > (72^ the 
higher loop orders are not suppressed relative to the one- 
loop contribution, and the whole approximation breaks 
down. 

In the electroweak theory the ratio A/52 is essentially 
given by the ratio of the Higgs and W boson masses. The 
Higgs boson has not been found and its mass is therefore 
unknown, but there is an experimental limit that its mass 
has to be at least 115 GeV. This means that X > g^, and 
the perturbative calculation breaks down. This is also 
the case in Type II superconductors. In both cases, one 
needs non-perturbative methods to study the properties 
of the phase transition. 



NON-PERTURBATIVE RESULTS 

In the global theory, one can be certain that no matter 
what non-perturbative effects there may be near the tran- 
sition, the transition must exist, because one can measure 
the value of (</>) = (f>a. It vanishes in the symmetric phase 
and is non-zero in the broken phase, and a function that 
does that cannot be analytic. In other words, (0) acts as 
an order parameter. 

When one considers non-perturbative effects in the 
gauge theory, one makes the simple but rather striking 
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observation that {(j)) cannot be an order parameter, be- 
cause it is not physically observable 0- This is because 
it transforms non-trivially under gauge transformations, 

W^)>-e"(-)(^(a:)). (28) 

If one averages over all gauge-equivalent configurations, 
the expectation value of any such quantity vanishes. 

The most reliable way of studying phase transitions 
non-perturbatively is to use lattice Monte Carlo simu- 
lations 0, and they show that the theory in Eq. H18|l 
still has a transition. At small the transition is 

of first order as our perturbative calculation indicated, 
but when A/g^ is larger, it becomes a continuous second- 
order transition. Even though one cannot use (0) as an 
order parameter, one can see which phase the system is 
in by studying the properties of the magnetic field. It is 
massless in the symmetric phase and massive in the bro- 
ken phase The latter is equivalent to the Meissner 
effect in superconductors. 

In non-Abelian theories, there are generally no mass- 
less fields, because the theories are confining. Therefore 
there is no qualitative difference between the "symmet- 
ric" and "broken" phases of the electroweak theory. Care- 
ful lattice Monte Carlo simulations have shown that when 
the Higgs mass is above 80 GeV, the two phases are actu- 
ally smoothly connected and there is no phase transition 
between them 9]. Because the true Higgs mass has to 
be in this range, it seems likely that there was never a 
true electroweak phase transition in the early universe. 
On the other hand, this conclusion may still change if 
the Higgs sector of the electroweak theory is more com- 
plicated than in the simplest version of the theory. 

Similarly, the GUT phase transition may have been a 
smooth crossover from one phase to another, provided 
that X/g^ was large. This is perhaps surprising, be- 
cause the broken phase of any GUT is supposed to be the 
SU(3)xSU(2)xU(l) phase of the Standard Model, and it 
contains a massless photon. However, the photon mass 
is not exactly zero in GUTs It is finite, but sup- 

pressed by exp(— Tgut/T"), and therefore the photon is 
massless for all practical purposes. Nevertheless, an arbi- 
trarily small mass is enough to make a smooth crossover 
possible. 

CONCLUSIONS 

We have seen that symmetries that are broken in the 
vacuum are typically restored at high temperatures, and 
that this leads to phase transitions. In global theories, 
the transitions are generally of second order, and in gauge 
theories they are of first order if the scalar self coupling 
A is small relative to the gauge coupling g^. However, in 
the opposite limit, the transition may become a smooth 
crossover between two phases, analogously to the transi- 
tion between water and vapour at very high pressures. 



At first sight, it seems there have been at least three 
phase transitions related to particle physics in the early 
universe: the deconfinement transition of QCD, the elec- 
troweak transition and the GUT transition. On closer 
scrutiny it turned out that it is likely that they may have 
all been merely smooth crossovers rather than true phase 
transitions in the mathematically rigorous sense. Never- 
theless, they may still have important consequences es- 
pecially in the form of defect formation. 



Lecture II: 
Defect Formation 

INTRODUCTION 

Topological defects are non-linear, time- independent 
solutions of field theories, which are made stable by the 
topology of the theory. They have been studied exten- 
sively in cosmology condensed matter physics |l2j 
and, at a more theoretical level, in the context of quan- 
tum or classical field theory. 

In cosmology, magnetic monopoles formed at the GUT 
phase transition would have disastrous effects later on in 
the evolution of the universe, and part of the motivation 
behind the inflationary theory jl^ was to wipe them out 
before that. Cosmic strings, on the other hand, are fully 
compatible with what is known about cosmology. If they 
were formed at some early stage, we may well be able to 
detect them with astronomical observations. 

Topological defects also exist in many different con- 
densed matter systems. The ones that are closest to de- 
fects of cosmological importance are vortices in superflu- 
ids and Abrikosov fiux lines in superconductors. 

It was pointed out by Kibble in 1976 yjj that topo- 
logical defects are generally formed at symmetry break- 
ing phase transitions and that we should therefore ex- 
pect that they were formed in the early universe, too. 
However, as discussed in Ref. 0, the known cosmolog- 
ical phase transitions are not simple symmetry breaking 
transitions, but they involve a breakdown of a local gauge 
symmetry. Nevertheless, as we shall see, the same con- 
clusion still holds: Topological defects are formed. 

In recent years, defect formation has been studied ex- 
tensively in condensed matter s^tems such as liquid crys- 
tals [illiSlli^uperfiuids [H US HI and super- 
conductors [23 Hi m ill. The first two of these in- 
volve global symmetries, and indeed, they seem to con- 
firm Kibble's theory. Superconductors have a U(l) gauge 
symmetry and would therefore be closer to the cosmolog- 
ical phase transitions. The superconductor experiments 
carried out so far have not produced conclusive results, 
but hopefully this will change soon. 
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In this lecture, I will first explain how topological de- 
fects form in transitions associated with a breakdown 
of a global symmetry. This is known as the Kibble- 
Zurek mechanism. Then, I will discuss gauge symmetries 
and show that the mechanism is actually rather differ- 
ent in that case. For a more comprehensive review, see 
Refs. mil. 

DISCRETE GLOBAL SYMMETRIES 

Consider first a real scalar field in a potential with a 
Z2 symmetry, 

- im202 + 1 (29) 

If is positive, the only minimum is at </> = 0, and the 
state of the system is symmetries under Z2. However, 
for rn^ < 0, there are two degenerate minima at (f) — 

If the field is initially exactly zero, it remains zero even 
if is negative. However, in real systems there are al- 
ways both quantum and thermal fiuctuations. Although 
the arguments we will use in this section are probably 
also valid for pure quantum fiuctuations, we assume that 
the temperature is non-zero and the dominant fiuctua- 
tions are thermal. Therefore, we will not consider any 
quantum effects. 

If either the value of <j) or its velocity is initially non- 
zero, it determines which minimum the system ends 
up in. If the initial condition is inhomogeneous, i.e., 
position-dependent, different parts of the system may 
end up in different vacua. On the whole, each minimum 
is equally probable, and therefore the system ends up in 
a configuration of positive and negative regions, which 
both cover half of the whole system. The boundaries 
between these domains are topological defects known as 
domain walls. The system is in vacuum on both sides of 
the wall, but because of continuity, the field must leave 
the vacuum inside the wall to interpolate between the 
positive and negative values. This means that the wall 
has positive energy per unit area. 

The most obvious question about the domain configu- 
ration is whether it has a characteristic length scale. In 
general, the initial fluctuations have a finite correlation 
length ^, which is defined by 

{mm) - {m)){m)) - e-i-~-^i/«. (so) 

This means that if the points x and y are further away 
from each other than the initial conditions at these 
points are independent. If we then assume that the do- 
mains form in a short enough time that no signal can 
travel from x toy before they have formed, there is a 50% 
chance that the system will end up in different vacua at 
these two points, in which case there would have to be (at 



least) one domain wall between them. Thus, one expects 
the typical distance between the walls to be roughly ^. 

Let us now assume that is initially positive but 
decreases with time. Initially, the system is in thermal 
equilibrium in the symmetric phase, but when reaches 
a critical value, a symmetry breaking phase transition 
takes place. Because the leading effect of a non-zero tem- 
perature is to change the effective mass term , we can 
equally well assume that is constant, but T decreases 
due to, say, the expansion of the universe or a contact 
with a colder heat bath. 

It is a general property of second-order phase transi- 
tions that the correlation length diverges at the transition 
point. This divergence can be characterised by a critical 
exponent ly, which is universal in the sense that it only 
depends on the universality class and not on any micro- 
scopic details, 

?(r)-^o(^^) (31) 

Typical values of v are around 2/3. 

If we could cool the system adiabatically so that it 
stays in equilibrium, the correlation length ^ would be 
infinite at the transition temperature T^. Consequently, 
the whole system could end up in the same vacuum and 
there would be no domain walls. 

However, this is not possible in practice, because of 
critical slowing down |29|. The relaxation time r, i.e., 
the time it takes for a small perturbation to equilibrate, 
also diverges exponentially, with exponent /i, 

r(r)-ro(^^^) ^ (32) 

This means that if we are cooling the system with a con- 
stant rate characterised by the "quench time" tq, 

T{t) =(^l-^^ Te, (33) 

the equilibrium correlation length and relaxation time 
diverge as 

e(0=eo(|^) \ r(0=ro(|^) \ (34) 

Once T(t) > \t\, any deviation from equilibrium will not 
have time to disappear before the transition. The fluctu- 
ations in the system at that time freeze out and survive 
all the way to the transition. In particular this means 
that the correlation length cannot grow significantly af- 
ter that time, because it is simply a property of these 
fluctuations. 

This means that the actual correlation length at the 
time of the transition is 

i = at), (35) 
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FIG. 5: The Kibble-Zurek mechanism is based on the ob- 
servation that the choice of the vacuum at points A, B and 
C separated by more than the correlation length ^ must be 
uncorrelated. 

where t is the freeze-out time determined by the equation 
T{t) = \t\. This is straightforward to solve and gives 

t^-iroT^y/^'+^\ (36) 
whereby the frozen-out correlation length is 

?" = ^o(^) ocra/(^+-). (37) 

As we saw before, this is the typical distance between 
domain walls formed in the transition. 

CONTINUOUS GLOBAL SYMMETRIES 

It is straightforward to generalise these arguments to 
systems with continuous symmetries. In the case of a 
complex scalar field with a U(l) symmetry, the set of 
possible vacua, i.e., the vacuum manifold, is a circle. The 
vacua are characterised by the phase angle (p of the field 
(j) — ve^'^. After the transition, the angle (p is uncorre- 
lated at distances longer than ^. 

Consider now three points A, B and C separated by 
more than ^ (see Fig. |5} . Each one has a random value 
of (fi, but field cj) is still continuous. It will therefore have 
to somehow interpolate between the three points, and if 
the points are not separated by much more than ^, it is 
most likely that the field follows the shortest path on the 
vacuum manifold, because that costs least energy. This 
assumption is known as the geodesic rule. However, there 
is a finite probability that if the geodesic rule is satisfied 
between A and B, between B and C and between C and 
A, the phase angle changes by 27r as one follows it from 
A to B to C to A. In that case, continuity tell us that 
the field has to leave the vacuum manifold and vanish 



at some point within the circle formed by these points. 
That point is a vortex. 

As with domain walls, the typical distance between 
vortices after the transition is ^, and their number density 
is therefore 

n^r^o,TQ^''^^'+^\ (38) 

It is straightforward to generalise these arguments to 
more complicated symmetry groups. In particular, one 
finds that the number density of magnetic monopoles is 

In principle, Eq. H38|) makes a prediction that can be 
tested in condensed matter experiments. However, the 
arguments that led to it neglected factors or order one in 
many places, and therefore even if one knew exactly the 
values of the parameters smd tq, the prediction could 
still easily be off by a factor of ten. A more robust pre- 
diction is the dependence on tq, but that would require 
carrying out experiments with a wide range of different 
cooling rates, and in many cases it is not possible to vary 
the cooling rate at all. 

An even more robust test of the Kibble-Zurek mecha- 
nism is given by the spatial distribution of vortices. Let 
us define Nyy{R) as the winding number around a circle 
of radius R, 

Nw{R) = 7^ </" df-^ip. (39) 

On average, this of course vanishes, but its variance 
(NwiR)^) is typically non-zero. After a phase transi- 
tion in a global theory, the circle passes through around 
domains that are uncorrelated with each other. 
Every time, the curve moves from one domain to an- 
other, the phase angle changes by some random amount 
A(/?. Labelling the domains by i, assuming the geodesic 
rule, and ignoring all factors of 27r, the total winding 
number is then 

Nw{R)~Y.^'f-- (40) 

i=l 

Because the changes A(^i are uncorrelated, the total vari- 
ance is 

= (A^2^i?/e«i?/e (41) 

More generally, one can define a scaling exponent vd by 

M 

{Nw{Rf) ocR^"". (42) 

The above result then implies that the Kibble-Zurek 
mechanism predicts vd = 1/4- This can be compared 
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with a fully random distribution, in which the sign Sa 
of each vortex a G {!,..., A^vort} inside the curve is in- 
dependent of all the others. In that case, the winding 
number is 



Nw{R) = 



(43) 



a=l 



and consequently 



(44) 



Thus, a random distribution of vortices would have vjj 
1/2. 



GAUGE SYMMETRIES 

In a gauge field theory, one cannot apply the Kibble- 
Zurck arguments as such, because the phase angle ip is 
not a gauge invariant quantity. Furthermore, while a 
global system tries to minimise phase gradients V(p in 
order to minimise energy, the condition for minimum en- 
ergy in a gauge theory is W(p -|- = 0. If the gauge 
field A is non zero, the phase if does not even attempt 
to become uniform. 

In the special case in which the magnetic field is zero 
everywhere, one can choose a gauge in which A — 0. 
Then, the gradients of tp cost energy, and one can expect 
the same arguments to hold as in a global theory. 

However, if one considers a non-zero temperature, as 
we have been doing, one cannot make this assumption, 
because it is not consistent with thermal fluctuations 31] . 
In fact, these fluctuations are nothing but thermal radi- 
ation, and in the symmetric phase, they can be approx- 
imated by the blackbody spectrum. To the extent that 
we can consider classical dynamics, we can therefore say 
that the initial conditions for the magnetic field are given 
by the Rayleigh- Jeans spectrum. A convenient way to 
express that is by defining the function G(fc), 



{2T:f5{k + k'). (45) 



The Rayleigh- Jeans distribution corresponds to G{k) = 
T. 

In the low-temperature phase, the long-wavelength 
part of the spectrum is suppressed by the Higgs effect. 



G{k) 



TP 



(46) 



where m-y is the photon mass. These same considerations 
apply to superconductors, as well, and in that case 
is the inverse penetration depth. 



Following Ref . |15| , one can calculate function 



of temperature, and to leading order one finds, 



p2 

— (T^ 

4 ^ 



T2) 



Assuming again linear cooling, one finds 



t 



(47) 



(48) 



We will now assume that as the system is cooled 
through the transition, the scalar field remains close to 
equilibrium. This is not always true, and if it is not, then 
one would expect something like the Kibble mechanism. 
However, here we are interested in the opposite limit in 
which the gauge field falls out of equilibrium first. We 
also assume that the mean field result Eq. H48I) is approx- 
imately valid. 

According to Eq. H4t)|l . the amplitude G{k) of a given 
mode decreases slowly at first, but when w fc, it starts 
to drop sharply toward zero. Using Eq. (|48|l . we see that 
this happens at time 



2Jf_ 

e2T2 



(49) 



After this time, the equilibrium value of G{k) decreases 
as 

2fc2 TQ 



G,^{k) 



e^Te t 



and we find 



rflnGeq(fc) 
dt 



(50) 



(51) 



This rate is faster at early times, and the fastest decay 
takes place at t « i^. Thus, 



dlnGeq(fc) 



dt 



1 

Tk 



2k^TQ • 



(52) 



This is the rate at which G{k) has to be able to decay for 
the mode to stay in equilibrium. We can see that long 
wavelengths (low k) have to decay faster. 

On the other hand, one would expect on physical 
grounds that the dynamics of longer wavelengths is slower 
and that they would therefore actually decay shower than 
short wavelengths. It is therefore inevitable that some 
very long-wavelength modes are unable to reach the nec- 
essary rate in Eq. (|52|l and fall out of equilibrium instead. 

We can be more specific if we assume that the dynam- 
ics is governed by the conductivity u and Ohm's law is 
valid, 



J = <yE. 

Maxwell's equations then imply 



(53) 



B 



-V X V X B + V X /= -V X V X S - o-B. (54) 
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In Fourier space (and ignoring vector indices), this be- 
comes simply 



B + (7B + k^B = 0, 



(55) 



which is nothing but a damped harmonic oscillator. If 
k < cr/2, it has an exponential solution 



B cx exp I 1 



(56) 



For G{k), this means that any deviation from the equi- 
librium value will decay no faster than at this rate, 



dlnG(fc) 
Jt 



> - 



2fc2 



(57) 



Thus, we find that the critical value kc, defined in such 
a way that modes with k < kc fall out of equilibrium, is 
given by 



1/4 



(58) 



One can then roughly say that as — > oo, modes with 
k > kc decay to zero, whereas those with k < kc re- 
tain their original amplitude. Thus, we would have the 
spectrum 



G(fc) 



k <^ kc^ 
0, k > kc- 



What this means is that the magnetic field has not dis- 
appeared completely. However, because the field cannot 
penetrate the broken, superconducting, phase, the mag- 
netic field must be confined into flux tubes, i.e., vortices. 
Thus, we have seen that in the gauge theory, vortices are 
formed by the non-equilibrium dynamics of the gauge 
field. This phenomenon is knows as flux trapping [27|, 
because these vortices originate in the thermal fluctua- 
tions of the magnetic fleld, which do not have time to 
decay and are instead trapped in vortices. 

In general, the winding number Nw and the magnetic 
flux <I> = / (PS ■ B are related by the flux quantum $0 = 
27r/e, and for the case of a circle of radius R in Eq. (jSHl, 
we have 



Nw{R) = «'(i?)/$o ~ e$(i?). 



(60) 



On average, this is zero, but as in the global case, the 
variance is non-zero and can be calculated from G{k), 

{Nw{Rf) « e^{HRf) 



d^xd^y{B,{x)B,{y)) 



"2 / d^xd^y I -^e^'= (^-^")G(fc)(61) 



d^k 
(2^ 




FIG. 6: A snapshot from a two-dimensional simulation by 
Stephens et al. [32i | shows clearly the clusters of positive 
(white) and negative (black) vortices formed in a gauge field 
theory transition. 



Using Eq. this is 



{Nw{R)' 



e^TcR, if i? > 1/kc, 
e^TcklR'^, if i? < l/kc. 



(59) The scaling exponent V£, is therefore 



1, at short distances, 
1/4, at long distances. 



(62) 



(63) 



Unfortunately, the long-distance value 1/4 is equal to 
what the Kibble- Zurek mechanism gives, and cannot 
therefore be used to distinguish between the two theo- 
ries. Because it depends on the dimensionality, the sit- 
uation is different in two-dimensional models that have 
been studied numerically j^ij 32] . Note, however, that 
a two-dimensional film in three-dimensional space, the 
typical setup of actual experiments, behaves like a three- 
dimensional system in this sense l33| . 

In three dimensions, the best way to distinguish be- 
tween the theories is the short-distance exponent v — 1. 
It means that the typical value of Nw{R) grows as i?^, 
i.e., it is proportional to the area of the circle. This is 
only possible if all the vortices inside the circle have the 
same sign, or in other words, if there are clusters of equal- 
sign vortices (see Fig.EJ. Eq. ll^ tells us that the radius 
of these clusters is « 1/fcc, and the number of vortices in 
each is ~ e^Tc/kc- Consequently, the number density 
per unit cross sectional area is 
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Here we have only discussed Abelian gauge field theo- 
ries, but similar arguments seem to apply to non- Abelian 
theories as well. In that case, one can also study the for- 
mation of magnetic monopoles js^l , which may shed more 
light on the monopole problem in cosmology. 

CONCLUSIONS 

We have seen that if a symmetry breaking phase tran- 
sition takes place in a finite time, topological defects are 
formed. In a global theory, they are produced by the 
Kibble-Zurek mechanism, which leads to a characteristic 
negative correlation between vortices that can be used to 
identify the mechanism. In a gauge theory, defects are 
also produced by thermal fluctuations of the magnetic 
field that are too slow to decay. The clearest prediction 
of this flux trapping scenario is the formation of clusters 
of equal-sign vortices. 

These theories are currently being tested in condensed 
matter experiments, in particular with superconductors. 
Although the calculations done in these paper refer to rel- 
ativistic field theories, the qualitative conclusions would 
be the same in condensed matter systems. 

I would like to thank the organisers for a very success- 
ful workshop, and the ESF COSLAB programme and 
Churchill College for financial support. 
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